Nonlinear Competition Between Small and Large Hexagonal Patterns 
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Recent experiments by Kudrolli, Pier and Gollub Q on 
surface waves, parametrically excited by two-frequency forc- 
ing, show a transition from a small hexagonal standing wave 
pattern to a triangular "superlattice" pattern. We show that 
generically the hexagons and the superlattice wave patterns 
bifurcate simultaneously from the flat surface state as the 
forcing amplitude is increased, and that the experimentally- 
observed transition can be described by considering a low- 
dimensional bifurcation problem. A number of predictions 
come out of this general analysis. 

PACS numbers: 47.20. Ky, 47.54. +r 



In recent years there has been considerable interest 
in pattern formation in spatially extended nonequilib- 
rium systems 0]. These studies have identified universal 
mechanisms for the formation of certain patterns. For 
instance, the origin of the ubiquitous hexagonal pattern 
can be traced to a symmetry-breaking instability of a 
spatially homogeneous state that occurs in the presence 
of certain Euclidean symmetries, which play a central 
role in determining the possible nonlinear evolution of 
the instability. In this letter we use similar ideas to in- 
vestigate a family of hexagonal and triangular patterns 
which arc born in the same symmetry-breaking instabil- 
ity that gives rise to simple hexagons, but have structure 
on two disparate length scales [[^jj]. Kudrolli et al. Jl]] 
call such structures "superlattice patterns" ; Fig. [I] gives 
two examples. Current interest in superlattice patterns 
is sparked by recent observations of their formation in 
experiments on parametrically excited surface waves |l]] 
and in nonlinear optics experiments as well as in a 
study of Turing patterns in reaction-diffusion systems || . 

In this letter we show that superlattice patterns can 
arise directly from the spatially homogeneous state via a 
transcritical bifurcation. Moreover, this occurs in generic 
pattern-forming systems for which there is only one un- 
stable wavenumbcr. We show that the triangular su- 
perlattice pattern differs from simple patterns such as 
hexagons and stripes, and complex quasi-patterns, be- 
cause it is characterized by both an amplitude and a 
phase which depend on the distance from the bifurca- 
tion point. In order to investigate the phase associated 
with this state we find that it is necessary to include 
high order resonant interaction terms in the amplitude 
equations. When these terms are neglected the prob- 



lem becomes degenerate and superlattice patterns such 
as those in Fig. [j] become just two isolated examples in a 
continuum of states that have varying degrees of symme- 
try. Thus high order terms are essential to understanding 
the formation of superlattice patterns. 






FIG. 1. On the left is a blow up of the experimental su- 
perlattice standing wave pattern flll (courtesy of Kudrolli, 
Pier and Gollub). We show a region with area ~ 1/30 of 
the circular container's cross-sectional area. Note that the 
pattern is periodic on a hexagonal lattice, and that it has 
triangular symmetry. On the right is a hexagonal superlat- 
tice Turing pattern obtained from a numerical integration of 
a reaction-diffusion system ^ (courtesy of Judd). 

While much of our analysis is quite general, our pre- 
sentation will focus on the recent experiments of Q], in 
which surface waves are parametrically excited by sub- 
jecting a fluid layer to a periodic vertical acceleration 
that involves two rationally-related frequencies. In the 
classic Faraday problem, which has single frequency forc- 
ing, the first modes to lose stability with increased accel- 
eration are subharmonic with respect to the vibration 
frequency [7 . With two frequencies, harmonic response 
can occur [3 j|. In particular, the triangular superlattice 
pattern in Fig. [I] was obtained with an acceleration of 
the form f(t) — a(cos(x) cos(6wi) + sin(x) cos(7wt + <j))). 
In the experiments, a hexagonal standing wave pattern 
is produced at the onset of instability, then, as the ac- 
celeration is increased, there is a transition to the su- 
perlattice pattern depicted in Fig. [j] Qj. Comparison of 
the spatial Fourier transform of the onset hexagons and 
the superlattice patterns reveals that the new state is 
formed from the nonlinear interaction of twelve promi- 
nent Fourier modes whose wavenumbers have the same 
modulus k c and lie on a hexagonal lattice with funda- 
mental wavenumber k c j\fl . In this letter we address the 
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transition between hexagons and superlattice patterns 
by considering a degenerate bifurcation problem akin to 
that which successfully explains the transition between 
hexagons and rolls in a variety of systems. A number of 
concrete predictions about the form, origin and stability 
properties of the superlattice patterns come out of this 
general symmetry-based analysis. 




FIG. 2. Intersection of the critical circle |k| = k c with the 
k— space hexagonal lattice generated by ki, k2. In this ex- 
ample a — 3, P = 2 in (^). The wave vectors ±Ki, . . . , ±K;j 
lie at the vertices of a hexagon, as do ±K4, . . . , ±Ke. 

The experimentally observed patterns are put into a 
theoretical framework by restricting to solutions that tile 
the plane in a hexagonal fashion. The time-periodic forc- 
ing in the Faraday problem leads us to a formulation in 
terms of a stroboscopic map. Since we seek spatially pe- 
riodic solutions, we express all fields in terms of double 
Fourier series; for example, the height of the free surface, 
at time t — mT, is 

C m (x) = Re( »n(mT) e '(«k 1+ n 2 k 2 ).x^ ) (1) 

n£Z 2 

where T is the period of the forcing. We obtain patterns 
periodic on a hexagonal lattice when the primitive vec- 
tors ki,k 2 satisfy |ki| = |k 2 | = k and ki • k 2 = — \k 2 . 
A standard ansatz in pattern selection problems is to set 
k = k c , where k c is the critical wavenumber of the insta- 
bility. However, here we want to investigate competition 
between small and large hexagonal patterns, which we do 
following the approach in jlj. Specifically, we choose k 
so that there are twelve wave vectors K n = niki + n 2 k2 
in (|l|) that satisfy |K n | = k c , i.e., so that there exists a 
co-prime integer pair n = (a, (3), a > (3 > 0, such that 

|aki + /3k 2 | = ky/a 2 + f3 2 - a(3 = k c . (2) 

Fig. H presents an example associated with a = 3, /? = 2, 
k c /k = \fl. Note that the lattice points on the critical 
circle lie at the vertices of two hexagons rotated by an 
angle 9 relative to each other. The angle, determined by 
the dot product of Ki = aki + /3k 2 and K4 = aki + {a — 

/3)k 2 , satisfies cos(6>) = Za^-aB+a*) ■ ^ on ^ tnose modes 



associated with one of the hexagons are excited, then 
the periodicity of the pattern is dictated by k c and one 
recovers the standard six-dimensional hexagonal lattice 
bifurcation problem. However, if all twelve modes on the 
critical circle are excited, then the period of the pattern 
is greater by a factor of k c /k. In order to simplify our 
presentation we now set a — 3,(3 — 2; the results for 
general a, (3 are similar. 

In formulating the bifurcation problem, we assume 
that the flat fluid surface loses stability to harmonic 
waves of wavenumber k c as a bifurcation parameter A 
is increased through zero. Thus the Faraday instabil- 
ity sets in when a Floquet multiplier fi crosses the unit 
circle at fi = 1. In our formulation there are twelve 
(real) Fourier modes in (|]) that are neutrally stable at 
A = 0; all others are damped. In this case the strobo- 
scopic map can be reduced, near the onset of instabil- 
ity, to a twelve-dimensional center manifold spanned by 
the critical Fourier modes. Let Zj(m), j — 1, ... ,6, be 
the complex amplitude, at t — mT, of the Fourier mode 
e jKj X , with Zj being the amplitude of e _lKj ' x . Here the 
Kj are labeled as in Fig. 0. 

The possible nonlinear terms in the stroboscopic map 
z(m + 1) = f (z(m)), z = (zi, . . . , zq), are restricted by 
the symmetries of the problem. The first component of 
f has the general form 

/1 = hi (u, q, q)zi + h 2 (u, q, q)z 2 z 3 

+ Cl z 2 z 2 z 6 + c 2 z 1 z 2 z 4 z 2 + 0(\z\ e ), (3) 

where u = (\z 1 \ 2 ,\z 2 \ 2 ,\z 3 \ 2 ,\z i \ 2 ,\z 5 \ 2 ,\z e \ 2 ), and q = 
(ziz 2 Z3, z^z^zq). The discrete hexagonal symmetries 
place further restrictions on the functions h\,h 2 and 
also determine the other components of f from fx [Q. 
At this point we note that if the resonant interaction 
terms, with coefficients C\ and c 2 , are neglected, then 
the phases 4>j associated with each Zj = Rje l ^ j only en- 
ter the problem through the values of the total phases 
$1 = <f>\ + <j) 2 + 4>3 and $ 2 = (j)4 + fa + cf)Q. As a result the 
cubic truncation of (^|) admits solutions in the form of two 
hexagons rotated relative to each other by « 22° (for 
a = 3, (3 = 2), and translated relative to each other by 
arbitrary amounts. Kudrolli, et al. |l| noted that only a 
restricted set of relative translations of the two hexagons 
leads to a superlattice pattern which retains triangular 
symmetry. Their phcnomcnological description did not, 
however, give a mechanism for the selection of a partic- 
ular shift. We now investigate an alternative description 
of the triangular superlattice pattern which has the ad- 
vantage that is in fact consistent with the inclusion of the 
high order resonant interaction terms. 

We focus on patterns that have at least three-fold ro- 
tational symmetry by restricting the bifurcation problem 
to the subspace z(m) = (u m , u m , u m , v m , v m , v m ), where 
u m and v m are complex Fourier amplitudes measured at 
time t = mT. First we recall some results about period- 
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one simple hexagons which satisfy (u m ,v m ) = (Re l,p ,0) 
or (0, i?e* v ). The amplitude R and phase (p obey 

= A + eRcos{3ip) + AR 2 + ■■■ 

= sin(3</?)(-e + 5i? 2 + •••), (4) 

where e, A and B are nonlinear coefficients that arise 
from the Taylor expansions of hi and h 2 in (||). The solu- 
tions of (|i|) that bifurcate from the origin at A = satisfy 
sin(3(p) = 0; the hexagons with <p = 0, are related by 
translations, as are the hexagons with (f = it, These 
two sets of hexagons, 'up- hexagons' and 'down- 

hexagons' (H~), form the two branches associated with 
a transcritical bifurcation. In non-Boussinesq convection, 
these states correspond to ones in which fluid is rising or 
falling at the centers of convection cells Gfl] . 




FIG. 3. Examples of patterns, periodic on a hexagonal 
lattice, that bifurcate transcritically from the flat state at 
A = 0. These are plots of appropriate superpositions of criti- 
cal Fourier modes with wave vectors ±Ki , . . . , ±Ks of Fig. ^[ 
The plots are, left to right, down-hexagons (H~) and super- 
lattice down-triangles (T~) (cf. Fig.Q). The critical wavenum- 
ber k c dictates the size of the small scale structure evident in 
the superlattice pattern. 

Next we consider patterns that involve all six complex 
Fourier modes zi, . . . , zq, with Zj — pe 1 ^ . We find that 
the amplitude p and phase ip satisfy equations of the 
form, cf. (|), 

= X + epcos{3i})+Ap 2 + ■■■ (5) 
= sin(3V0 (-e + Bp 2 + ■ • •) + (c x - c 2 )p 3 sm(20) + • • • 

In this case there are two distinct types of solutions that 
bifurcate from the origin at A = 0; those with phase 
ip = Q, 7T and those with ip w ±7r/3, ±27r/3. Note that the 
phase ip associated with the latter solutions depends on 
the amplitude p, which in turn depends on the distance 
from the bifurcation point. The solutions with tp = 0, tt 
have hexagonal symmetry; in p| they are referred to 
as "super hexagons" . As with simple hexagons, super 
hexagons bifurcate transcritically with the two parts of 
the branch satisfying ip = and ip = tt, respectively. 
The solutions satisfying ip ±7r/3, ±27r/3 have only tri- 
angular symmetry, and are new. The triangular solutions 



with ip w 27r/3 and ip w — 7r/3 form the two parts of a 
transcritical branch; a rotation of these patterns by tt 
changes the sign of ip. These solutions have structure 
identical to the superlattice harmonic state described in 
[jlj : compare the triangular superlattice pattern in Fig. |l| 
with that in Fig. ||. Moreover, we know that the solutions 
with cos(3(p) > in ([!]) and cos(3i/>) > in (||) bifurcate 
in the same direction from the origin, as do the branches 
with cos(3</j), cos(3-0) < 0. 

Finally, we address the experimentally observed tran- 
sition between the patterns in Fig. ^. Since all of the 
known primary solution branches are unstable at bifurca- 
tion, due to the presence of the quadratic nonlinear term 
in (||), we assume that the coefficient e of the quadratic 
term satisfies |e| <C 1. This allows us to investigate stable 
small amplitude solutions, and the transitions between 
them. If we truncate (0) at cubic order, thereby neglect- 
ing the high-order resonant terms with coefficients c\ , C2, 
then the superlattice patterns are at best neutrally sta- 
ble. Specifically, while we expect a (multiplicity two) unit 
multiplier associated with translations of the pattern, the 
extra multiplier (also of multiplicity two) results from a 
symmetry of the truncated equations that allows a rela- 
tive translation of the two rotated hexagons that make 
up the pattern. As noted above, this symmetry is broken 
when the resonant terms are included; the Floquet multi- 
plier /i then moves off the unit circle. In particular, p, > I 
for the super hexagons if 4ci -I- 5c2 < Q . Analogously, 
we can show that p < 1 for the triangular superlattice 
pattern if 4ci +5c2 < 0. Thus if the superlattice patterns 
are neutrally stable when the high-order resonant terms 
are neglected, then we expect one of them is stable and 
the other is in fact unstable. This is consistent with the 
experimental observations of Kudrolli, et al. [3] who only 
observe the triangular superlattice pattern. 

In Fig. ^, we present part of a bifurcation diagram 
associated with the stroboscopic map. Specifically, we 
keep all terms through cubic and the essential quintic 
terms: 

fl = (l + \)zi + ez 2 z 3 + (ai|zi| 2 + a 2 \z 2 \ 2 + a 2 \z 3 \ 2 )z! 
+ (a 4 |z 4 | 2 + a 5 |z 5 | 2 + a 6 \z e \ 2 )zi (6) 
+c\z\z\z§ + c 2 ziz 2 Z4zl , 

and assume that the cubic coefficients satisfy the follow- 
ing two inequalities: 

a 1 + 2a 2 < — \a 4 + a 5 + a 6 | (7) 

Hal - a 5 ) 2 4- (a 4 - a 6 ) 2 + (a 5 - a 6 ) 2 
a 2 - «i > y , 

with 4ci + 5c2 < 0. Conditions (0) ensure that superlat- 
tice patterns eventually supersede the simple hexagons, 
and that the stripe pattern is never stabilized. 
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FIG. 4. Part of the bifurcation diagram associated with 
the bifurcation problem (^) for — 1 <C e < 0. Bifurcation 
points are indicated by solid circles. Only branches of simple 
hexagons (H ) and superlattice triangles (T ) are indicated; 
solid lines correspond to stable solutions. Branches produced 
in secondary bifurcations, and unstable primary branches are 
not shown. 

This example demonstrates that it is possible to re- 
produce the type of transition reported in within the 
framework of a generic bifurcation problem on a hexago- 
nal lattice. Moreover, a number of predictions come out 
of this general analysis. Specifically, we find that both 
the transition between the the flat state and the small 
hexagons, and the transition between the hexagons and 
the superlattice pattern are hysteretic, with the region 
of bistability between the hexagons and the superlattice 
pattern exceeding that between the hexagons and the 
flat state. While there is no clear detection of hysteresis 
in the experiments, the transition from hexagons to the 
superlattice pattern is reported to occur via domains of 
the superlattice pattern growing, with increased accelera- 
tion, until they fill the container Jl| . Coexisting domains 
of different patterns, possibly with slowly moving fronts 
between them, is a common manifestation of bistability 
of patterns in extended systems. Another prediction that 
comes out of our general analysis is that for e < the 
transition is between the down-states depicted in Fig. [|, 
while if e > the transition is between up-states. While 
in Rayleigh-Benard convection there is a clear distinc- 
tion between up- and down-states, this distinction may 
be a subtle one in the Faraday experiment since stan- 
dard imaging techniques rely on reflection of light from 
the surface of the fluid, so that surface peaks are not dis- 
tinguishable from troughs. The nonlinearity inherent in 
the imaging may also present challenges to extracting the 
amplitude-dependent phase that we predict is associated 
with the triangular superlattice patterns. 

While the foregoing analysis shows how it is possible 
to achieve the experimentally observed transitions within 
a general framework, many of the specific predictions of 
the analysis still need to be tested. Furthermore, it would 



be very interesting to calculate the nonlinear coefficients 
in (^) from the hydrodynamic equations to determine 
whether the minimal inequalities (0) are satisfied for the 
experimental parameters. Such computations are quite 
involved in the case of two-frequency forcing of viscous 
fluids. However, much progress has been made in the 
case of very low viscosity fluids by Zhang and Vinals 
p"T[ , who employ a quasi-potential approximation to the 
Navier-Stokes equations to simplify the calculations. 

Finally, we note that all our calculations have assumed 
that the solutions are strictly periodic on a hexagonal 
lattice. This is certainly an appropriate model for obser- 
vations of superlattice patterns, but as noted in jjj quasi- 
patterns are observed for other experimental parameters. 
A mechanism for favoring quasi-patterns by a resonant 
interaction between two bifurcating states with different 
horizontal wavenumbers was proposed by Edwards and 
Fauve ||, and investigated for a Swift-Hohenberg type 
of model by Lifshitz and Petrich Whether a simi- 

lar resonance mechanism is responsible for the spatially- 
periodic superlattice patterns is an intriguing, open ques- 
tion. 
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